Orbital Liquid in Perovskite Transition-Metal Oxides 
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We study the effects of the degeneracy of the e g orbitals as well as the double exchange interaction 
with t2 g spins in perovskite transition-metal oxides. In addition to the spin field Si, the isospin field 
Ti is introduced to describe the orbital degrees of freedom. The isospin is the quantum dynamical 
variable, and is represented by the boson with a constraint. The dispersion of this boson is flat along 
(n/a, n/a, k z ) (a: lattice constant) and the other two equivalent directions. This enables the orbital 
disordered phase down to low temperatures. We interpret some of the anomalous experiments, i.e., 
optical absorption and d.c. resistivity, in the low temperature ferromagnetic phase of Lai-zSr^MnC^ 
with x > 0.2 in terms of this orbital liquid picture. 
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The colossal negative magneto-resistance observed in 
transition-metal oxides with perovskite structure, e.g., 
the double exchange system La 1 _ 2; Sr 2 ,Mn03, has revived 
the interest in these systems The valency of Mn 

ion is Mn 3+ for x = 0, whose electron configuration is 
(i2g) 3 (e £( ) 1 , and all the spins are aligned ferromagneti- 
cally due to the strong Hund coupling. Because of the 
strong on-site repulsion, the double occupancy of e g or- 
bitals is forbidden and the system is a Mott insu- 
lator. When one La is replaced by Sr, one hole is intro- 
duced to Mn and Mn 3+ turns into Mn 4+ . These doped 
holes, which contribute to the conduction, can not be 
described in terms of the one-body theory. For a single 
band case, which is relevant to high-Tc cuprates, exten- 
sive studies have been focused on these doped holes to a 
Mott insulator. One of the promising approaches is res- 
onating valence bond (RVB) theory |7|J§], where the spin 
and charge are separated. In the orbital-degenerate case, 
this idea can be generalized also to the orbital degrees of 
freedom as shown below B. 

Experimentally the orbital ordering has been estab- 
lished in the low hole-concentration region, i.e., x ~ 0.0 
, where the system is insulating with the A- type an- 
tiferromagnetic long range ordering (AFLRO) As 
x increases the system becomes more and more conduc- 
tive, and finally shows the metallic conduction below the 
ferromagnetic transition temperature T c . In this metallic 
state (x > 0.2), there are several anomalous features. 

(1) The optical conductivity shows a narrow coherent 
peak up to around 0.02eV followed by a broader "Drudc 
like" band up to around leV Jl2| |. The integrated oscil- 
lator strength for this two-component Drude absorption 



changes down to the very low temperature where the 
ferromagnetic moment already saturates, which suggests 
that other degrees of freedom still remain active. 

(2) The photoemission spectra shows only a small dis- 
continuity at the Fermi energy Ep followed by a gap like 
behavior |l3|-|l5f| . On the other hand, the specific heat is 
not enhanced with 7 ~ 5mJ/K 2 mol Jig] . 

(3) The anisotropy of the conduction and the spin exci- 
tation, which is expected with the orbital ordering, is not 
observed even at low temperatures |T(^Jl7|] . 

(4) There is no symptom of the Jahn- Teller (JT) dis- 
tortion in neutron scattering experiments fllTI . The dis- 
placements of the oxygen ions are independent of the 
temperature across the ferromagnetic transition temper- 
ature. 

It is noted that the on-site repulsive interaction within 
each orbital plays no role in the ferromagnetic state be- 
cause of the Pauli's exclusion principle. Then it is sug- 
gested from the above features that the orbital degrees 
of freedom play some roles in the low temperature phase. 

In this paper we study the fluctuation in the orbital 
degrees of freedom in perovskite transition-metal ox- 
ides. These degrees of freedom are represented by the 
isospin T, which is the quantum dynamical variable with 
T = 1/2 for the e g orbitals. We found that its fluctu- 
ation has a low dimensional character. Based on this, 
we propose that the isospin is still disordered, i.e., liquid 
state, in the ferromagnetic state of Lai-^Sr^MnOa with 
x > 0.2. We have done some numerical simulations based 
on this idea, and their results are at least encouraging. 

In order to study fluctuation in the orbital degrees of 
freedom, we adopt the following Hamiltonian. 



H = e d J2 <a V + E (*« '4y^ a + h.c) + H e ^ e + -K £ d\ 1<T (a) aa A^ ' ^ ■ (1) 

Here it is assumed that the 2p orbitals of oxygen have been integrated over, and only the d orbitals of transition- 
metal ions are considered. The operator d\ creates an electron with spin cr(=f , |) in the orbital 7(= a, b) at site 

i. The transfer intensity tJJ between 7 orbital in site i and 7' orbital in the nearest-neighbor site j is calculated by 
considering the overlap integral between the d and p orbitals. We choose the up and down isospin state in the orbital 
space as a = (3z 2 — r 2 ) and b — (x 2 — y 2 ) orbitals, respectively. Then tj 7 is explicitly written by 
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for fj = fi ± x, fj = fi ± y and r) = fi ± z, respectively. The constant to is positive, and depends on the distance 
between the d and p orbitals. The electron-electron interaction iJ e _ e is explicitly written by 

H e - e = UJ2 n nT^7l + U 'Y1 n ™ nih + 1 d \ao d \bo' d l ao' d l bo > ( 3 ) 

27 i i,(T,(7' 



where n i7CT = 4 7(T ^ 70 - and n. n = X) ff n» 70 -. Tne fourth- 
term in Eq. (|l|) is the Hund coupling between e g and 

t2 g spins and S^ 2s is the t2 g spin with S = 3/2. The 
electron-electron interaction H e - e in Eq. (||) is rewritten 
by using the spin operator S = |E-, ff( /4'^^'V' 
and the isospin operator T = i E 77 '<t ^o-f^^'^V for 
the orbital degrees of freedom. First notice the identities 
S 2 = |n-f ^ 7 n 7T n 7i +2^ a -5 6 ,andf 2 = \n-2§ a -§ h - 
n a rib + ~ ^ 7 n 7 |n 7 |, where n — n a + rib. Then ff e _ e is 



rewritten as H e _ e = Sf + a®T? + PS ia - S ib } , 

where aW = \U + \U' - \l, a® = W - ±1, and 
(3 = — I (U — U' — I). Because (i) j3 is order of / and 
is much smaller than and and (ii) the large U 1 
forbids the simultaneous occupancy of a and b orbitals, 
we will neglect the Si a ■ Sib term. Since a'*' and a'*) 
are both positive, the electron-electron interaction tends 
to induce the spin and isospin moments. By introducing 
two kinds of Stratonovich-Hubbard auxiality fields, the 
partition function is given by 

D^\r)D^\r)y-I drL , (4) 



with 

L = £ di~fa(d T + e d )d lia + £ ifij d h^ d j7'^ + h - c -) + K ^2$i ' ^ + L t 2g 

i,7,c (ij) ,(7,7,7' * 

-2£(a« ft>.3 t +aW • 3 ) + £(>VT 2 + 2 ) , (5) 



where L t2g is the Berry phase term of the ti g spins. Eq. 
(||) describes the e g electrons moving in the background 
of two fluctuating fields 2a {s) 4> {s) and 2a {t) $ {t) . We in- 
troduce at this stage the approximation which is appro- 
priate for the strong correlated case. In this case the 
magnitude of the local fields a^\<p^\ and \(p^ are 
much larger than the transfer intensity to, and the elec- 
tron is forced to be aligned in the directions of and 
at each site. Thus, it is convenient to rotate the spin 
and isospin axes in each site, in order that its z-axis co- 
incides with </)( s ) and 4>^\ This is accomplished by the 
unitary matrices U^' and for the spin and orbital 
spaces, respectively, which transform the fermion opera- 
tor as / <7 v/ = {U> a) ) a , a (UW) , di . In this rotating 
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frame the fields cf>^ and <ft^ are pointing in the direction 
of +z, and accordingly the density of states for /-fermion 
are divided into four bands separated by the energy gaps. 
When the concentration of the e g electrons is one per each 
transition-metal ion, only the lowest band, which corre- 
sponds to er =f and 7 = a, is occupied and the system 
becomes a Mott insulator. We keep only this lowest band, 
because the holes are doped to it and only this is impor- 
tant when the low energy excitations are concerned. We 
introduce the spinless and orbital-less hole operator hi as 
hi = // a | . The virtual transition processes to the higher 
bands cause the exchange interaction terms Lj, i.e., so 
called J term in the t-J type models |is|| . Then the ef- 
fective Lagrangian up to the second-order with respect 
to the electron transfer is obtained as follows, 

, a + (U^d T Ut%^j 

+ J2Hu^dM s) )u -£ + Lj , (6) 

» (ij) 



where the chemical potential /x/, is determined by the 
condition (h\h t ) = x. The term 3(U^d T U^ s) ) II comes 
from Lt 2g in Eq. (Q). The Berry phase terms of the 



original electron di ia generates those for the rotated 
fermions hi, spins, and isospins. Here we introduce the 
spinor boson zf^ 1 =* [zffl, zty] to represent the unitary 
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where $ s) /l$ S) | 



and J2cr \ z ia\ 2 = 1- Correspondingly 



for spins are written as (U^ d T U^)^ 



J2af3 Z ld K a )a(i Z i(i 

we introduce 2} for C/^* . Then the Berry phase terms 

= T z {s) *d z (s) 
and its coefficient should be regarded as the spin quan- 
tum number 2S. Then in the undoped case (x = 0), 
the quantum numbers of the spin and isospin are S = 2 
and T = 1/2, respectively. Then it is expected that the 
quantum fluctuation is stronger for isospin than spin. 

Now the original electron operator di ia with the con- 
straint of no double occupancy is expressed as di ia — 

h\ z i-y z ia ) which is the generalization of the slave- fermion 
formalism to the orbital degenerate case. Then we call hi, 



zl', and z[ s J holon, isospinon, and spinon, respectively. 
In analogy with the spin-charge separation in high-Tc 
cuprates H, one can consider the orbital-charge separa- 
tion in this formalism, which we discuss below. Then the 
effective transfer intensity ty is given by 
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To study the nature of the orbital fluctuations in the 
ferromagnetic state, we derive the effective action for the 
isospinon by employing the mean-field approximation. 
We first consider the kinetic energy hiUjhj as- 

suming that xto 3> J where J is the typical exchange 
energy and is of the order of t 2 /a^ or t^/aS^. Then 
in the limit of strong correlation, we have the concen- 
tration region where J /to <C x <C 1, which we are now 
interested in. The exchange terms Lj will be discussed 
later. We replace the fermion operators by their average 
value, i.e., (hihi) — x, {hihj) — —x, and quench the mag- 
netic fluctuation by putting z\ s J =* [1,0]. The constraint 

J2-y \ z< hJ\ 2 = 1 is imposed on average by introducing the 
chemical potential —A. Then the effective Lagrangian of 
the orbital excitation is obtained as 
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The eigenvalues for each k are given by e± (fc) = A — 
xhf±(k) where f±(k) = -{c x + c y + c z ) ± [c 2 x + c 2 + 
c 2 — c x c y — c y c z — CzCx] 1 / 2 , with c x = cos(ak x ) etc. The 
minimum of this energy are given by the flat dispersion 
along the axis (it /a, it /a, k z ) and the other two equiva- 
lent directions. This situation is quite contrast to the 
spin case where the dispersion relation is proportional to 
cos(afc a; ) + cos(afcj / )-|-cos(afc 2 ). The electron in x 2 — y 2 or- 
bital can not hop perpendicular to the xy-plane because 
the overlap with the oxygen p orbital is vanishing due to 
the symmetry. ( Only (aa) component is non-zero in the 
third matrix in Eq. (^).) 

This flat dispersion leads to an important conse- 
quence. The chemical potential A is determined by 

(T,j z i^* z iy) = ^E 7 fc n s(4 t) (fc)) = 1 where n B is 
the bose distribution-function. Because of this flat dis- 
persion, the chemical potential —A is always negative 
at finite temperature although it becomes exponentially 
small below the effective bose condensation temperature 
T* ~ xto. Then the orbital long range ordering does not 
occur down to low temperatures. 

One may think that the physical quantities behave 
similarly to the orbital ordered state when the chemi- 
cal potential is very small and z^> is almost bose con- 
densed. However this is not the case because there are 
three-branches of the low lying fluctuations, i.e., along 
(jr/a,Tr/a,k z ) and two equivalent directions. In order to 
show it, we generate the random configuration of z^> 
at very low temperature where the chemical potential is 
already very small. In this case only the static com- 
ponent of z^' is important, which is obtained by the 
random numbers generated by the Gaussian function 



exp[— e± (fc)/T]. Then is obtained through Eq. (Q), 
and we diagonalize the holon Hamiltonian for that con- 
figuration of the transfer. We study the density of states 
and the optical conductivity cr^ui) of the holon system 
by averaging over 50 random samples for the cubic lattice 
of 8 x 8 x 8. The results for the optical conductivity is 
shown in Fig. I]. If one look at U/s, they fluctuate in sign 
and magnitude violently, which explains the incoherent 
nature of <7h(u>). 



FIG. 1. The optical conductivity of the holon (7h{u) cal- 
culated in the 8x8x8 system. The temperature is chosen 
to be T = O.lto- 

Now we consider the gauge- field fluctuations. There 
are two kinds of U(l) gauge field corresponding to z^ 
and z'*', respectively. The spatial components of which 
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are given by the phase of the order parameters, i.e., 
(*&N3> = Ik^^S^lexp^), (r = s,t A, A' = 
cr, 7). The time components, on the other hand, are the 
Lagrange multiplier fields which impose the constraint 
I z Ia' > I 2 = I n * ne ferromagnetic state, z^ is bose 
condensed, and becomes massive, while a'*) remains 
massless in the orbital disordered state. More impor- 
tantly the Ioffe-Larkin composition rule |Q applies also 
the orbital-charge separated state, and the optical con- 
ductivity <j{u>) is given by cr(u>) — ah{oj)at(uj) / (ah(u>) + 
a t {ijj)). Then the characteristic energy of the orbital fluc- 
tuation, which is much smaller than that of holons, can 
manifest itself in a(tu) through <Jt{uj). The photoemission 
spectrum, on the other hand, is given by the convolution 
of those for holon and isospinon, because the Green's 
function G in space-time representation is the product 
G(r, t) — Gh(r,r)Gt(r,T). According to the flat disper- 
sion of the isospinon, the holon dispersion is averaged 
over and the Fermi edge will be smeared out as observed 
experimentally. 

Here we discuss several possible effects which break the 
flat dispersion and/or make the orbitals to be long range 
ordered. One possibility is so-called J terms, i.e., Lj in 
Eq. (||). It is noted that the transfer matrix tij = (tjj ) 
does not commute with any of the Pauli matrices a a . 
Then there is no continuous rotational symmetry in the 
isospin space. Neglecting the terms of the order of xJ, 
Lj in the ferromagnetic phase has the same form as the 
first term in Eq. (0) with Za (%h ) being replaced by T z 
(T x ) and xt by -%/(U'-I) @. Note that only T z and 
T x appear in Lj. The dispersion has the same form as 

z^-boson, i.e., ej±(k) = — u t ,°_ I f±(k). Then the mini- 
mum of the energy are given by the flat dispersion along 
the axis (ir/a, it/a, k z ) and the other two equivalent di- 
rections. For the lowest dispersion along (jv/a,Tt/a,k z ), 
T z = and only T x is nonzero, which competes with 
the fact that the kinetic energy prefers nonzero T z . As 
before the ordering of the isospin is suppressed because 
there are three branches of the flat dispersion with the 
role of z's being replaced by T's. 

Applying the mean-field approximation to Lj, which 
gives the two-body interactions between the z&' -bosons, 
Lj does not generate a dispersion for the flat branch be- 
cause flatness comes from the fact that the hopping of the 
x 2 — y 2 electron along the z-axis is forbidden. We also 
note that the orbital angular momentum is quenched in 
e g orbitals, and the spin-orbit interaction is ineffective. 

The orbital liquid, on the other hand, should be 
very sensitive to the anisotropy between x-, y-, and z- 
axis. In this respect, in the layered materials, e.g., 
Lai_2,Sri +;r Mn04 and La2-2zSi'i + 2xMn207, the orbital 
liquid should be absent. Experimentally the incoherent 
band extending to u> ~ leV seems to be missing in cr(u>) 
in these materials |^o| . Another test on our scenario is the 
effects of the uniaxial pressure. We expect cr(u>) becomes 
anisotropic and becomes sharp in the more conductive 
plane. 
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